Experimental evidence for Tayler instability in a liquid metal column 
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In the current-driven, kink-type Tayler instability (TI) a sufficiently strong azimuthal magnetic 
field becomes unstable against non-axisymmetric perturbations. The TI has been discussed as a 
possible ingredient of the solar dynamo mechanism and a source of the helical structures in cosmic 
jets. It is also considered as a size limiting factor for liquid metal batteries. We report on a liquid 
metal TI experiment using a cylindrical column of the eutectic alloy GalnSn to which electrical 
currents of up to 8 kA are applied. We present results of external magnetic field measurements that 
indicate the occurrence of the TI in good agreement with numerical predictions. The interference 
of TI with the competing large scale convection, resulting from Joule heating, is also discussed. 



The last years have seen a number of liquid metal ex- 
periments on various magnetohydrodynamic instabilities 
with relevance to the origin and the action of cosmic mag- 
netic fields [l|. Dynamo action has been observed in the 
large scale liquid sodium experiments in Riga, Karlsruhe, 
and Cadarache The helical version of the magne- 
torotational instability (MRI) has been evidenced in the 
PROMISE experiment [2,], and further experiments in 
Maryland [3] and Princeton (Ef are devoted to the inves- 
tigation of the standard version of MRI. What is missing 
yet in the liquid metal lab is any evidence of the Tayler 
instability (TI) [6]. This is in remarkable contrast to 
the vast experience in plasma physics where the (com- 
pressible) counterpart of TI is better known as the kink 
instability in a z-pinch [7|, i.e. the limit of the Kruskal- 
Shafranov instability when the safety factor goes to zero. 
In astrophysics, TI has been discussed as a possible ingre- 
dient of an alternative, nonlinear stellar dynamo mecha- 
nism (Tayler- Spruit dynamo [8]), as a generation mech- 
anism for helicity [9] , and as a possible source of helical 
structures in galactic jets and outflows [10]. 

A particular motivation to study the TI in a liquid 
metal arises from the growing interest in large-scale liquid 
metal batteries as cheap means for the storage of highly 
intermittent renewable energies. Such a battery would 
consist of a self-assembling stratification of a heavy liq- 
uid half-metal (e.g. Bi, Sb) at the bottom, an appropri- 
ate molten salt as electrolyte in the middle, and a light 
alkaline or earth alkaline metal (e.g. Na, Mg) at the 
top. While small versions of this battery have already 
been tested [ll| , for larger versions the occurrence of TI 
could represent a serious problem for the integrity of the 
stratification. In a recent paper [12] we have proposed a 
simple trick to avoid the TI in liquid metal batteries by 
just returning the battery current through a bore in the 
middle. By the resulting change of the radial dependence 
of B^p ir) it is possible to prevent the condition for (ideal) 
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FIG. 1. Experimental set-up. Left: Scheme with liquid metal 
column and fluxgate sensors positioned along the vertical axis 
and the azimuth. Right: Photograph of the central part of 
the experiment. 



TI, d{rB^{r))/dr > 0. Without such a provision, TI is 
expected to set in at some finite electrical current in the 
order of kA. The precise value is a fucntion of various 
material parameters, since for viscous and resistive fluids 
TI is known [l3| to depend effectively on the Hartmann 
number Ha = B^R{cr / {puy^'^ ^ where R is the radius of 
the column, a the electrical conductivity, p the density, 
and u the kinematic viscosity of the fluid. 

In this paper, we present experimental results that con- 
firm the numerically determined growth rates of TI [l3| 
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as well as the corresponding prediction that the critical 
current increases monotonically with the radius of an in- 
ner cylinder [T2l.[l^. 
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FIG. 2. Magnetic fields measured along the vertical axis, for 
ri=12.5 mm and /=7 kA. (a) Time dependence of Bz at the 
11 fluxgate sensors, (b) Power spectral density in dependence 
on time and wavelength. Note the dominance of short wave- 
length signals in the initial phase, and the dominance of the 
long wavelength signal at later times, (c) Detailed time evo- 
lution of the PSD for four significant wavelengths, showing a 
transition to an exponential growth around t ^150 s, which 
is later stopped at t ^280 s. 



The central part of our TI experiment (Fig. 1) is an 
insulating cylinder of length 75 cm and inner diameter 
10 cm, filled with the eutectic alloy GalnSn which is liq- 
uid at room temperatures. The physical properties of 
GalnSn at 25°C are: density p = 6.36 x lO^kg/m^, kine- 
matic viscosity u = 3.40 x 10~^m^/s, electrical conduc- 
tivity a = 3.27xl0^(l^m)~^ At the top and bottom, the 
liquid metal column is contacted by two massive copper 



electrodes of height 9 cm which are connected by water 
cooled copper tubes to a DC power supply that is able to 
provide electrical currents of up to 8 kA. By intensively 
rubbing the GalnSn into the copper, we have provided 
a good wetting making the electrical contact as homoge- 
neous as possible. 

While in later experiments it is planned to use Ul- 
trasonic Doppler Velocimetry (UDV) in order to mea- 
sure directly the axial velocity component along the z- 
axis, for the first experiments we have decided not to use 
any inserts that could possibly disturb the homogeneous 
current going from the copper electrodes to the liquid. 
Therefore, for the identification of TI we exclusively rely 
on 14 external fluxgate sensors that measure the vertical 
component of the magnetic field. Eleven of these sen- 
sors are aligned along the vertical axis (with a spacing of 
7.5 cm), while the remaining three sensors are positioned 
along the azimuth in the upper part, approximately at 
15 cm from the top electrode. The distance of the sensors 
from the outer rim of the liquid metal column is 7.5 cm. 
This rather large value, which is certainly not ideal to 
identify small wavelength perturbations, has been chosen 
in order to prevent any saturation effects of the fluxgate 
sensors in the comparably strong azimuthal field of the 
axial current. 

The main goal of our experiment is to study the in- 
fluence of the electric current through the fluid on the 
growth rate and on the amplitude of the magnetic field 
perturbations. This is done without any insert, as well as 
for two different radii of an inner non-conducting cylin- 
der, ri=6 mm and 12.5 mm, for which we expect a mono- 
tonic increase of the critical current. 

For the particular case with ri=12.5mm and / = 7kA, 
Figs. 2 and 3 illustrate the measured Bz data in depen- 
dence on the vertical position and on the azimuth, re- 
spectively, and the procedure of data analysis. In both 
cases, the time stamp and the field values were set to 
zero a few seconds after the switch-on process of the cur- 
rent had been finalized. Figure 2a shows the subsequent 
evolution of Bz measured at the 11 vertical positions. 
Most of the data in Fig. 2a show a collective long-term 
trend, whose source might be increasing convection, but 
possibly also some geometric changes (e.g. by thermal 
expansion of constructional parts) which could influence 
the projection of the strongly dominant azimuthal field 
component on the measured z-component. 

Apart from this long term trend, the data of Fig. 2a 
contain much more details on the vertical dependence 
that can be extracted. For this purpose, we first subtract 
(at every time instant) the mean value and the linear 
trend in z-direction and compute then the Power Spec- 
tral Density (PSD) of the remaining signal. The resulting 
PSD, in dependence on time and wavelength, is shown in 
Fig. 2b. At the beginning we observe the simultaneous 
growth of three modes with short wavelengths, while later 
a long wavelength mode becomes dominant. The time 
dependencies of four exemplary modes with wavelengths 
18 cm, 24 cm, 36 cm, and 72 cm are shown separately in 
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FIG. 3. Magnetic field data measured along the azimuth, 
for n=12.5 mm and /=7 kA. (a) Time dependence of Bz at 
the 4 fiuxgate sensors. The m=l character is clearly visible, 
(b) Squared magnetic field amplitude of the m=l mode, (c) 
Angle of the m=l mode. 



Fig. 2c. Evidently, after some initial transient (lasting 
approximately until 150 s), the growth of all short wave- 
length modes acquires an exponential character, and we 
can read off the corresponding growth rates in the time 
interval between 200 s and 280 s (see the dotted lines in 
Fig. 2c). Note the similarity of the growth rates of the 
18 cm and 36 cm mode which indicates the merging of 
next-neighbor cells (of the pure TI) as it is known from 
numerical simulations of the combined action of TI and 
(thermal) convection due to internal heating. 

This regular exponential growth of the short wave- 
length modes stops suddenly at t=280s, when the long 
wavelength mode (here with wavelength 72 cm) becomes 
dominant. Most interestingly, the growth rate of this 
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FIG. 4. (a) Measured growth rates and their numerical pre- 
dictions for three different radii of the inner cylinder, in de- 
pendence on the current, (b) Measured squared magnetic field 
strength of the m=l mode, and their piecewise linear fits. The 
maximum squared fields for ri=0 mm, 6 mm, and 12.5 mm are 
61.0 /iT^ 543.3 /iT^ and 1668.1 /iT^ respectively. 



72 cm mode after the transition becomes quite similar 
to the growth rate of the 36 cm mode shortly before the 
transition. This may indicate a sudden doubling of the 
wavelength with increasing temperature and convection. 

Now we turn to the discussion of the azimuthal de- 
pendence of the induced ^^-perturbations, as measured 
by the four sensors around the cylinder. The behaviour 
shown in Fig. 3a clearly indicates a non-axisymmetric 
(m=l) mode of the field with a growing amplitude at 
the beginning and some rotation and further growth at 
later times. The time evolution of the squared amplitude 
(5^(0°)-5^(180°))2 + (5^(90°)-5^(270°))2 of this m=l 
mode is presented in Fig. 3b, and its corresponding angle 
in Fig. 3c. 

In Fig. 4 we compile the dependence of various quan- 
tities on the radius of the inner cylinder and on the cur- 
rent through the fluid. The growth rates of the 18 cm 
mode, as extracted from the periods with clear exponen- 
tial growth (see Fig. 2c for one example), are compared 
with the numerically predicted growth rates in Fig. 4a. 
Despite some scatter of the data, we observe a quite rea- 
sonable agreement with the numerical predictions. 

The second quantity of interest is the saturation level 
of the magnetic field. Similar to many other instabilities, 
for TI one would expect a more ore less linear behaviour 
of the saturated squared magnetic field above the critical 
current, according to B'^ ^ {I — /cdt)- However, the spe- 
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FIG. 5. Comparison of the determined critical currents with 
the numerical predictions for the m=l mode. The numerical 
prediction for the m=0 (sausage) mode is indicated partly. 



cific problem of our liquid metal experiment is that the 
saturation of the TI does not rely on an intrinsic back- 
reaction process, but on the extrinsic stabilizing effect of 
the increasing large scale convection. With this caveat 
in mind, in Fig. 4b we try to identify the critical cur- 
rent by plotting the squared magnetic field of the m=l 
mode at the inflection point of the corresponding curves 
(indicated, for this example, by the point in Fig. 3b). 

Starting at this inflection point we typically observe 
also a drift of the angle of the field pattern. We hy- 
pothesize that this "mode locking" starts when the long 
wavelength mode becomes dominant and aligns to some 
preferred azimuthal direction (determined, e.g., by slight 
geometric imperfections of the experimental set-up). The 
observed magnetic field level at this inflection point is 
thus yet determined by the TI alone. Insofar, it is rea- 
sonable to assume a monotonic, if not linear, dependence 
of the quadratic field strength (taken at the inflection 
point) on (/ — /crit) to characterize the TI. 

By applying a piecewise linear fit to the observed 
squared field intensities (see Fig. 4b), we can determine 



the critical current at the crossing point of the two 
straight lines. Figure 5 shows the critical currents de- 
termined that way and the error bars (with a 95 per cent 
confidence level) , in comparison with the numerically de- 
termined ones. Despite the fact that the inflection-point 
criterion for the determination of the saturated mag- 
netic field strength is certainly debatable, we obtain again 
a reasonable correspondence with the predicted critical 
currents. 

The investigation of a current driven instability in a 
column of a liquid metal has revealed a rather complex 
behaviour, comprising (after some initial transient) an 
exponential growth of the anticipated TI mode followed 
by a large-scale convection dominated by Joule heating. 
From the sharp bends of the amplitude evolution of the 
short wavelength modes (see Fig. 2c) it becomes apparent 
that the growing intensity of convection limits the further 
growth of the TI. This is a drawback of the present exper- 
imental configuration, since it prevents any investigation 
of the intrinsic saturation mechanism of TI. The theo- 
retical picture for this mechanism relies on the on-set of 
strong turbulence (due to the TI) leading to a radially de- 
pendent turbulent resistivity (/3-effect, see [15]) which, in 
turn, would result in a modified radial current distribu- 
tion that is just marginally stable against TI. The inves- 
tigation of this interesting saturation mechanism would 
require a significant weakening of the role of Joule heat- 
ing, either by choosing a horizontal orientation of the 
column (which, however, introduces an explicit symme- 
try breaking), or by an increase of the radius, and/or by 
the replacement of GalnSn with liquid sodium, for which 
TI is expected to occur already around 1 kA. Such a large 
scale sodium experiment, in which TI and MRI can be 
studied together, is planned for the future. 
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